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Abstract. It is commonly accepted that a neutron star is produced, when a massive 
star exhausts its nuclear fuel and ends its life in a core-collapse supernova explosion. 
This scenario is confirmed by the detection of pulsars, which are believed to be rapidly 
spinning neutron stars, in the central regions of many supernova remnants. Neutron 
stars and their associated supernova remnants should therefore have the same ages. As 
expected, the age of the Crab Pulsar, the first to be connected with a supernova remnant 
(the Crab Nebula), can be inferred from its current spin period and its derivative and 
indeed has about the same age of the supernova remnant, that was produced from a 
historically recorded supernova explosion in 1054. However most neutron stars appear 
to be much older than the ages of their associated supernova remnants, a puzzle not 
yet understood. Another puzzle is that so far no convincing evidence has been found 
in favor of magnetic field decay in neutron stars, that is predicted in most models of 
neutron stars. Here we show convincing evidence of magnetic field decay in some 
young neutron stars, and that the magnetic field decay can alter their spinning behaviors 
significantly such that these neutron stars appear much older than they really are. 



1. Magnetic Field Decay of Neutron Stars 

A rapidly spinning neutron star loses its rotational energy by magnetic dipole radiation 
(Pacini 1968) and outflow of charged particles accelerated in its magnetosphere (Gol- 
dreich 1969), i.e., to - -Kur' l{2n)^ , where a» is its angular frequency and K oc II, 
in which M is its dipole magnetic moment and / is its moment of inertia; its surface 
dipole magnetic field strength can thus be estimated as, B = 3.3 x \0^^{PP)^l^ , where 
P and P are its spin period and period derivative, respectively. Then assuming ^ = 0, 
the neutron star's age can be calculated as, T^, - {P - Po)/2P, or if P » Pq, Tg = P/2P, 
commonly referred to as its characteristic or spin-down age. More generally, the spin- 
down law can be written as, d> = -Koif{2n)^~'^ , or, P - KP^~". Defining a measurable 
quantity, the so-called braking index, n\^ = 2 - PP/P^, it is easy to show (Lyne 2006), 
«b = « - kP^~" IP. Assuming n - 2>, n\, < 2> indicates K > 0, ox vice versa (Lyne 
1975, Chanmugam 1995, Lyne 2004, Lyne 2006, Chen 2006). However the number of 
neutron stars with reliably measured n\^ is currently very small (Chen 2006). 

Many studies on the possible magnetic field decay of neutron stars have been done 
previously, based on the observed statistics of P and P of neutron stars. Some of these 
previous studies used the spin-down ages as indicators of the true ages of neutron stars 
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and found evidence for their dipole magnetic field decay (Pacini 1969, Ostriker 1969, 
Gunn 1970); however as we will show in the following, their spin-down ages are nor- 
mally significantly different from their true ages. Alternatively other studies relied 
on population synthesis that includes complicated observational selection effects, un- 
known initial parameters of neutron stars, and not completely understood models of 
radio emission of pulsars (Holt 1970, Bhattacharya 1992, Han 1997, Regimbau 2001, 
Gonthier 2002, Guseinov 2004, Aguilera 2008, Popov 2010). Consequently so far there 
is no consensus on if and how neutron star's magnetic field decays (see (Harding 2006a, 
Ridley 2010, Lorimer 2011) for reviews). 

The true age of a neutron star can be given from the age of the supernova remnant 
physically associated with it, since they both were produced during the same supernova 
explosion. A supernova remnant expands and interacts with its surrounding interstel- 
lar medium, and thus its age, Tsnr, may be calculated by modeling its morphology at 
the current epoch. The ages of neutron stars obtained this way are independent of the 
properties of neutron stars and thus unbiased in principle. In figure 1 (upper panel) we 
show the distribution of T'snr/T's; clearly for many of them Tsnr/Ts 1, i.e., many 
neutron stars look much older than their associated supernova remnants. The distri- 
bution is possibly bi-modal, divided by Tsnr/^s ~ 0.5 ; we will return to this point 
later. In figure 1 (lower panel) we show the correlation between Tsnr/T^s and the dipole 
magnetic field strength B. The strong and significant positive correlation suggests that 
low values of B at the current epoch is the major cause for Tsm./Ts 1. This im- 
mediately suggests that magnetic field decay plays a significant role in making neutron 
stars appear older than they really are, at least for neutron stars younger than about 
one million years. This general idea was suggested previously (Lyne 1975, Geppert 
1999, Ruderman 2005), but so far has not been compared between data and models 
statistically. 

With the true ages of neutron stars available, we can then study the correlations 
between their ages and their other two observables normally available with each pulsar, 
namely, P and P, as shown in figure 2. Assuming heuristically that all neutron stars 
were bom with the same properties, the trends shown here can be explained as the 
evolutionary tracks of neutron stars. These trends suggest that, as neutron stars become 
older, their periods become longer, their period derivatives becomes smaller, and their 
dipole magnetic field becomes weaker. One may argue that the observed trend of their 
dipole magnetic fields may be a mathematical artifact, because B and P are positively 
correlated mathematically. However this argument does not hold, because B and P 
are also positively correlated in the same way mathematically. It is actually physically 
more meaningful to conclude that it is the magnetic field decay of neutron stars that 
causes P decreases rapidly with time. 

Several physical mechanisms have been proposed for magnetic field decay in 
neutron stars, e.g., ohmic dissipation, Hall effect, and ambipolar diffusion. We can 
therefore have (Goldreich 1992, Heyl 1998), ^ = -S(-^ + + t-), in which 



of 10^^ G, denotes the core temperature in units of 10^ K, p^uc = 2.8 X 10^^ g cm ^ 
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Figure 1 . Upper Panel: Distribution of the ratio between the age of a supernova 
remnant (Tsnr) and the spin-down age (T^) of the associated neutron star. The verti- 
cal dashed and dotted Hnes correspond to Tsnr - Ts and Tsnr - 0.5 Ts, respectively. 
Lower Panel: Correlation between the ratio and the dipole magnetic field strength. 
The solid, dotted, and dashed lines correspond to linear fitting results (in double 
logarithm scales) with all data points, the left-most point removed, and the five left- 
most points removed, respectively. In each case a strong and significant positive 
correlation is found. 
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Figure 2. Correlations between the true ages (Tsnr) of neutron stars and their 
periods (P), period derivatives (P), and surface dipole magnetic field strengths (B = 
3.3 X IQ^'^iPpy^^). The straight lines are linear fitting results (in double logarithm 
scales). 
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is the nuclear density, and L5 is a characteristic length scale of the flux loops through 
the outer core in units of 10^ cm. 

Assuming does not vary or only varies slowly with time (Pons 2009), the ohmic 
dissipation produces exponential decay of magnetic field, B = Bqc"'^'^'^, where tj is 
referred to as the decay time scale. The other two mechanisms all produce power-law 
decay, B = So(f/;o) where the decay index a = 1 or a = 0.5 for Hall effect or 
ambipolar diffusion (for T «c 10^^ years), respectively. For the exponential decay, the 
true age of a neutron star is given by (Lyne 1975) = 5Tdlog(l + ^^)- For the 

power-law decay, the true age of a neutron star is given by Tp = ^^^^^^Tg (for a < 0.5 

and Tp » to), or T^logiT^/to) = j^T, (for a = 0.5), or Tp - 'i£2zllT,(Tp/to)^"-' 
(for or > 0.5), respectively (in the above to and Bo are the age and dipole magnetic field 
strength of a neutron star when it started its current phase of magnetic field decay). 
Because of the coupling between to and a when a > 0.5, we will first discuss the case 
a < 0.5 in the following. We stress that only P - KP^~" and n = are assumed in 
deriving the above relations. 

In figure 3, we compare predictions of the exponential or power-law magnetic field 
decay with the observed correlation between and Tsnr of neutron stars; for conve- 
nience we take n = 3 for the moment and we will discuss more general cases later. 
Clearly the exponential decay model cannot describe the data points for most young 
neutron stars. The power-law decay model can describe the data points if different pop- 
ulations of neutron stars are allowed to take different values of a. The dipole magnetic 
field strength of each neutron star is shown in the same panel, suggesting that neutron 
stars with weaker dipole magnetic field tend to have larger a; this is a generic prediction 
of the model if their initial dipole magnetic fields have similar strengths. 

In figure 4, we show the distribution of a, derived from figure 3 by counting the 
number of neutron stars in each interval of a. A striking bi-modal distribution is ob- 
vious, indicating that about half of neutron stars have no or very weak magnetic field 
decay (a ~ 0) and the other half have very rapid magnetic field decay (a « 0.5). 
This bi-modal distribution is actually another manifestation of the bi-modal distribu- 
tion shown in figure 1 (upper panel); however here the two peaks are well separated, 
indicating that the power-law description of their dipole magnetic field decay is appro- 
priate and physical. If we take a = 0.5, then Tp \og{Tp/to) - T^, i.e., it is possible to 
reproduce all observed data points with Tp/T^ > 1 by adjusting to. For a > 0.5 and 
a - 0.5 X 0, (Tp/to)^"'^ ~ 1 and thus a w 0.5 can also reproduce all observed data 
points with Tp/T^ > 1. Therefore we conclude that the ambipolar diffusion (a - 0.5) 
dominates the magnetic field decay for nearly half of the young neutron stars associated 
with supernova remnants, implying that fambip < ^ohmic and fambip < '^haii- Clearly the 
first condition can be easily met for physically plausible parameters of neutron stars. 
The second condition, i.e., fambip < ^haii requires Sq > 6 X 10^^ G, about the dipole 
magnetic field strengths of these neutron stars without significant magnetic field decay, 
i.e., TsNR « Tg. 

It is notable that there is one neutron star (PSR B 1757-24) with < Tsnr- In 
this case, A" > is required, taking n = 3 as above (please do not confuse n with nb 
here and throughout this paper). This agrees with the suggestion that its non-recovered 
glitches have significantly increased its dipole magnetic field (Lin 2004). In addition, 
the distribution peak at Ts^^ - T^, shown in figure 1 (upper panel), requires n = 3, if 
it = for these neutron stars. Therefore n = 3 assumed above is reasonable. On the 
other hand, the general expression of the spin-down age is „ = j^P/P- Because P 
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Figure 3. Observed correlation between the spin-down Ts and true (Tsnr) ages of 
neutron stars, compared with the predictions of the exponential or power-law mag- 
netic field decay models. The solid line corresponds to Tsnr = ^s. The dot-dashed 
lines correspond to Tg - irj log(l + 2^), with - 10"*, 10^, and 10^ (in top-down 
order), respectively. The dashed lines correspond to Tp = (1 - 2a)Ts (for a < 0.5 
and T-p s> fo), with a - 0.499, 0.40 and 0.10 (in top-down order), respectively. The 
dividing values of a are chosen in such a way that the numbers of neutron stars in 
the three regions (in three different colors) are roughly the same, for visual clarity 
only. The strength of the dipole magnetic field of each neutron star is also illustrated 
with the radius of the circle proportional to log(B/(10'°(G))). 
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Figure 4. The histogram shows the distribution of a, derived from figure 3 by 
counting the number of neutron stars in each interval of a. The dashed curve plots 
Tsnr/Ts = (1 - 2a). The horizontal dotted line corresponds to Tgm. = O.ST^. A 
bi-modal distribution is obvious, indicating that about half of neutron stars have no 
or very weak decay (a x 0) and the other half have very rapid decay (a ~ 0.5). This 
bi-modal distribution is actually another manifestation of the bi-modal distribution 
shown in figure 1 (upper panel). 

clearly increases with time statistically, as shown in figure 2 (upper panel), P « Pq is not 
plausible. Therefore n » 3 is required for these neutron stars with Tsnr <s: if there 
is no significant magnetic field decay; this is clearly also not plausible. Alternatively, 
n < 3 makes the spin-down age (Pj „) of a neutron star even longer than assuming n - 3; 
in this case more or even all neutron stars have Ts „ > T^snr, i e., A' < for more or 
even all of these neutron stars. We therefore conclude that for any reasonable values of 
n, at least some of these neutron stars must have experienced significant decay of their 
dipole magnetic fields. 



2. Discussion 

It has been believed that magnetic field decay powers the persistent X-ray emission of 
magnetars (Thompson 1995, Thompson 1996). It has also been suggested that their 
magnetic field decay plays a significant role in their thermal evolution (Pons 2009), 
and thus has significant implications to understanding the neutron star's equation of 
state. Here by comparing the spin-down and true ages of young neutron stars, we 
find convincing evidence for magnetic field decay in about half of them. For magnetic 
field decay dominated by the ambipolar diffusion (a = 0.5), it is easy to show that 
Tg « 2.6 X 10"^ TT^^'^- Taking the typical values of S12 - 1, ? ~ 5 x 10"^ yr, and L5 ~ 1, 

its core temperature is predicted to be roughly ~ xlO^ K, corresponding to a surface 
temperature (Arras 2004) of ~ 6 x 10^ K and consistent with observations (Pons 2009). 

It however remains puzzling why these neutron stars discussed here, that are ob- 
served as normal radio pulsars, have experienced such diverse evolutionary routes: the 
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magnetic fields increase in very few neutron stars, remain nearly constant in about half 
of them, or decrease rapidly in others. Perhaps these neutron stars were bom with dif- 
ferent physical properties (Han 1997); understanding the initial conditions of neutron 
stars will inevitably advance our understanding of stellar evolution, supernova explo- 
sion, and the equation of state of neutron stars significantly. 
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